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Abstract 

We analyze mathematical and physical properties of a previously introduced [J. 

Phys. A47, 115302 (2014)] family of 17(4) coherent states (CS). They constitute 
a matrix version of standard spin U{ 2) CS when we add an extra (pseudospin) 
dichotomous degree of freedom: layer, sublattice, two-well, nucleon, etc. Appli¬ 
cations to bilayer quantum Hall systems at fractions of filling factor v = 2 are 
discussed, where Haldane’s sphere picture is generalized to a Grassmannian pic¬ 
ture. We also extend Wehrl’s definition of entropy from Glauber to Grassmannian 
CS and state a conjecture on the entropy lower bound. 
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1 Introduction 

Spin-s, Bloch, SU( 2) or atomic coherent states (CS) were introduced by Radcliffe [lj, 
Gilmore mm and Perelomov [I]. They, together with standard (Glauber) CS [5], ac¬ 
curately describe the physical properties of many macroscopic quantum systems like in: 
quantum optics, Bose-Einstcin condensates (BEC), two-level systems, superconductors, 
superfluids, quantum Hall effects, etc; see e.g. [El El El El HQ1 El] for a selected panorama of 
applications (not only in quantum mechanics, but also in engineering). In particular, the 
ground state of many physical systems undergoing a quantum phase transition (QPT) is 
well described by a CS. Actually, it was Gilmore who introduced an algorithm [12] . which 
makes use of CS as variational states to approximate the ground state energy, to study 
the classical, thermodynamic or meanfield limit of some algebraic quantum models. This 
algorithm has proved to be specially suitable to analyze the phase diagram of Hamiltonian 
models undergoing a QPT like: the Dicke model of atom-field interactions P31 M, ESI US], 
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Bose-Einstein condensates m, the Lipkin-Meshkov-Glick model UBl US EDI EQ, vibron 
model for molecules [22l [23l [24] . etc. In some quantum phases, the ground state is in fact 
a (parity) symmetry adapted CS or “Schrodinger cat” [22JI2BJ- 

The generalization from U( 2) to U(N ) CS is quite straightforward for the fully sym¬ 
metric (canonical) case, related to the complex projective CP^ -1 = U(N)/[U(N — 1) x 
U( 1)] quotient space. However, general Grassmannian cosets Gj/ = U(N)/[U(N — k) x 
U(k)], related to other less symmetric Young tableau arrangements, are a bit more in¬ 
volved and perhaps less known or explicitly worked out. Actually, we are interested in 
the noncanonical chain of subgroups 

17(4) Z> 1/(2) x U (2) D U{ 1) x U{ 1). 


The corresponding matrix representation is useful, for instance, when adding a pseudospin 
(layer, two-well, nucleon, sublattice, etc) to the spin degree of freedom. For example, in 
bilayer spin (namely, quantum Hall) systems, pseudospin (pspin for short) is introduced 
in order to emphasize the spin SU( 2) symmetry in the, let us say, bottom (pspin —1/2) 
and top (pspin 1/2) layers |27j, ;28] ESj; here pspin rotates when particles are transfered 
from one layer to the other. Also, as noted long time ago by Heisenberg and Wigner 
[3QJ, U( 4) (which assigns isotopic spin —1/2 to protons and 1/2 to neutrons) is an ap¬ 
proximate symmetry of nuclear forces when neglecting electromagnetic forces. Moreover, 
pspin in graphene makes reference to the two triangular sublattices of the two-dimensional 
honeycomb structure. 

We shall make use of the bilayer quantum Hall (BLQH) picture all along this paper to 
provide a physical interpretation of some of the mathematical structures that will arise. 
In fact, spin-s CS already appear in Haldane’s sphere picture [31j for the (monolayer) 
fractional quantum Hall effect, in which 2s makes reference to total magnetic flux through 
the surface in units of the Dirac magnetic flux quantum <f> = 2t\ A/e; the spin s is also 
related to the “monopole strength” in the sphere § 2 . This image has to do with the 
“composite fermion” picture of Jain [32], which explains the existence of fractional values 
for the filling factor v (number of electrons per magnetic flux quantum penetrating the 
sample). The generalization of Haldane’s sphere § 2 = U{2)/U{1) 2 from monolayer to 
bilayer quantum Hall systems is the Grassmannian Gg = t/(4)/Z/(2) 2 (we simply write 
Go, since no confusion will arise) for fractional values of the Elling factor v = 2/A, with 
the isospin A related to the number of magnetic flux quanta piercing each of the two 
electrons (also the monopole strength in G 2 ); see [281 l29j for more information. 

The organization of the paper is as follows. Firstly in Section [2] we briefly review the 
simpler case of U( 2) CS, which allows us to better understand the more involved U( 4) 
case. In Sec. [3] we introduce different bases of U (4) operators, the carrier Hilbert space 
of the unirrep associated to the phase space f/(4)/G(2) 2 and the corresponding set of 
CS. We show that Grassmannian CS can be seen as a matrix generalization of Bloch CS 
with a generalized binomial distribution (Sec. 3.1). In Sec. 3.2 we compute operator 


symbols and relate them to order parameters in the BLQH jargon, studying the so called 
spin, pspin and canted quantum phases [27]. Using the CS (Bargmann) representation 
(see [33] for a review), in Sec. 3.3 we also compute the Husimi function Q^ of any state 
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0 and extend Wehrl’s definition of entropy [3T] from Glauber to Grassmannian CS. We 
conjecture a lower bound for this extension of Wehrl’s entropy (and an upper bound for 
the Husimi’s second moment) which is attained for any Grassmannian CS. Finally, Section 
Q] is left for conclusions. 


2 Bloch U(2)/U(l) 2 spin-s coherent states: a brief 


Let us consider the carrier space T~L S of a ( 2 s + l)-dimensional irrep of SU{ 2 ) spanned by 
the orthonormal basis vectors {|s, k), k = —s ,..., s}, which are eigenstates of the angular 
momentum third component J 3 with eigenvalue k. From now on we shall remove the spin 
label s from the basis vectors and simply denote | s,k) = \k). Spin-s CS can be obtained as 
an exponential action of angular momentum ladder operators J± — J\ ± iJ 2 on (namely) 
the lowest-weight state \k) — \ — s) as 


e zJ + 


■s) _ El=-sMz)\k) 

(1 + 

\z\ 

l 2 ) s (1 + 1 

\z\ 

2)s 


( 1 ) 


1 /2 

with <Pk(z) — ( s+ \.) 7 z s+k a set of monomials in the complex variable z verifying the 
closure relation (the Bergmann kernel K 2s ) 


y, <Pk{^)<Pk{z) = (1 + z'z) 2s = K 2 s (z', z ). (2) 

k=—s 


The complex label z is also usually parametrized as z — tan(6 l /2)e 1 ^, related to the 
stereographic projection of a point ( 6 ,0) (polar and azimuthal angles) of the Riemann 
sphere § 2 = U (2) jU (l ) 2 onto the complex plane. The probability of measuring the angular 
momentum third component k in the CS \z) follows a binomial distribution 


\(k\z )\ 2 


\Mz )\ 2 

(1 + \z\ 2 ) 2s 



(3) 


with parameter cos 2 (0/2). This binomial structure is also apparent from another familiar 
Fock-space representation [equivalent to 0 ] of spin-s CS as a two-mode, a and b, Bose- 
Einstein condensate (BEC) of iV = 2s particles, given by the expression (|0 f) denotes the 
Fock vacuum) 




1 / + za t 

vW v^rry 


2s 


\ 0 F ). 


(4) 


In bilayer quantum Hall (BLQH) systems -spin-frozen case- and in Bose-Einstein conden¬ 
sates in a double well potential, the operators cd and b 1 create (flux) quanta in layers/wells 
a and b, respectively {29i, 17] . In these physical contexts, the polar angle 9 is related to the 
population imbalance between wells/layers (see [29] and later on this section), represented 
by the operator J 3 = |(cba — tfb) in the usual Jordan-Schwinger boson realization for spin 
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J„ = 1 (a* 6t) (Tjj (a by, where denote the usual three Pauli matrices (/i = 1,2,3) 
plus the 2x2 identity matrix cro, and 2J 0 = cda-t-h^h is the total number of particles. The 
azimuthal angle (p is the relative phase between the two spatially separated wells/layers, 
an important magnitude to codify qubits in these systems. 

Spin-s CS Q are normalized, but not orthogonal in general, as can be seen from the 
CS overlap 




K 2s (z',z) 
K s (z', z')K s (z, z) ’ 


(5) 


written in terms of the Bergmann kernel 0. Spin-s CS constitute an overcomplete 
set and fulfills the resolution of unity 1 = J §2 \z)(z\dp(z, z), with measure dp(z,z) = 
= sin 9d6d(j) (the solid angle element). 

We shall make extensive use of the CS (Fock-Bargmann) representation 'h(^) (the 
“symbol”) of any state \ip) G T-L s , defined as 


'h(z) = K s (z,z){z\ip). 


( 6 ) 


For example, from 0 and 0, the basis states | ip) = | k) are represented by the monomials 
(Pk(z) °f degree s+k, whereas general states \ip) = Yhk C A^) are represented by polynomials 
T(^) = °f degree < 2s in the variable z. The relation 0 makes the space 

"H s (S 2 ) of polynomials of degree 2s in the variable 0 a reproducing kernel Hilbert space 
[that is, a Hilbert space of functions T in which pointwise evaluation T (z) is a continuous 
linear functional] with measure dp, s (z,z) = dn(z,z)/K 2 s(z,z). Inside this CS picture, 
angular momentum operators Jj are represented by differential operators J t , namely 

J'i = Z— - s, J+ = -Z 2 — + 2 sz, J- = — , (7) 


so that the following identity holds: J^{z) = K s (z,z){z\Ji\'ip). This relation is useful for 
technical calculations; for example, taking \ip) = \z'), we can easily compute CS matrix 
elements as 


(z\Ji\z') 


J l K 2s {z,z') 

K s (z,z)K s (z',z')' 


( 8 ) 


In particular, the diagonal elements (Jj) = (z\Ji\z) (also called “lower” or “contravariant” 
symbols of Jj) are simply (Jj) = K 2 f{z,z)J’iK 2 s{z,z). In the same way, CS expectation 
values of operator products can be computed as (JjJ*) = K^^z, z)Jj{JiK 2s (z, z)). For 
example, in the context of BEC/BLQH systems commented after 0 , the quantities 


(Js) 


z] 2 -l 
z I 2 + 1 


—s cos 9, (J 2 ) 


(Js ) 2 + 


2 s |^| 2 

(|^| 2 + 1) 2 


s 2 cos 2 9 + 



(9) 


are related to population imbalance (J3) and its fluctuations \J (Jf) — ( J3) 2 between 
wells/layers. As commented at the introduction, these CS expectation values are im¬ 
portant quantities to compute energy surfaces and to study the classical, mean-held or 
thermodynamical limit of spin systems undergoing a quantum phase transition (QPT). 
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Actually, spin CS saturate the Heisenberg uncertainty relation (J/)(Jf) — K^) 2 an d 
are therefore “close to classical”. Semiclassical properties of quantum spin states given 
by a density matrix p are also described by the CS expectation value Q(z,z ) = (z\p\z), 
also called Husimi quasiprobability distribution, lower symbol or Q-function of p. For 
pure states, p = |-0)(0|, we have Q{z,z) — l(^|'0)| 2 = | v h(^)| 2 /Jl 2 s (z,z). If (0|0) = 
1, then Q is normalized according to J g2 Q(z, z)dp(z, z) = 1. Wehrl’s entropy W = 
f s2 Q(z, z) In Q(z, z)dp(z, z) measures the area occupied by a state -0 in phase space, 
which is proved to be minimal, W m \ n = 1 — l/(2s + 1), when -0 is itself a CS |-0) = | z') 
(see [35| for a recent proof, which was conjectured in 1978 by Lieb [36]). Equivalently, 
the second moment (“inverse participation ratio”) M = f §2 Q 2 (z, z)dp(z, z) is maximum, 
M max = 1/2 + 1/(2 +8s), when 0 is a CS. In the large spin limit [37] [L] we have W m \ n —* 1 
and M max —^ 1/2, thus recovering the extremal values attained by ordinary Heisenberg- 
Weyl (harmonic oscillator) coherent states, as proved by Wehrl [33]. We shall also con¬ 
jecture extremal values of W and M for f/(4)/t/(2) 2 CS in the next section. 

To conclude this section, let us remark that all these information measures (Renyi- 
Wehrl entropies) have been of important use to characterize QPTs of several molecular 
and optical models [38] [23] [39] [2U 120] 00] and topological insulators [H] 02] 03] . 

There are many more interesting mathematical properties and physical applications 
of U(2)/U(l) 2 CS, but we just want to give here a brief account enough to introduce and 
compare with the more involved C(4)/£/(2) 2 CS in the next section. 


3 Grassmannian £/(4)/£/(2) 2 isospin- A coherent states 

We shall add an extra pseudospin (pspin) degree of freedom to spin, which will lead to 
some matrix generalizations of the standard spin-.s CS expressions of the previous section. 
Looking for physical applications, pspin may refer to a “layer, well, sublattice, nucleon, 
etc” extra degree of freedom that comes into play. 

To emphasize the spin-pspin (combined into “isospin”, to use the same nomenclature 
as in [27]) structure, we denote the U( 4) generators in the fundamental representation 
by the sixteen 4x4 matrices r= <r pspm <g) c+ pm , p, v — 0,1, 2, 3 (we shall use Greek 
letters /i, v for four-indices and latin letters i , k, l = 1, 2, 3 for the “spatial” part). General 
U( 4) unirreps arise in the Clebsch-Gordan decomposition of a tensor product of N four¬ 
dimensional (fundamental, elementary) representations of 17(4), for example: 


N=2\ 


□ ( 


.□ = 



2A 











© 





N= 2A 


or [1] <g> • • • <g> [1] = [N] © [A, A] ® ... (10) 


The best known case is the fully symmetric unirrep [first term in the direct sum decom¬ 
position (10)] corresponding to the Young Tableau of shape [N] (one row of N boxes) 
with dimension (N + 1 )(N + 2 )(N + 3)/6. Here we are interested in those unirreps cor¬ 
responding to the Young Tableau of shape [A, A] with two rows of A = N/2 (for even N) 
boxes each [second term in the direct sum decomposition (10)]. The dimension of the 
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(a) v = 2 


(b) v = 2/3 



Figure 1: Composite particles (bosons or fermions), seen as “fat electrons” bound to A 
(even or odd) magnetic flux lines, explain fractional quantum Hall effect at filling factors 
v = n/X. We represent n — 2 electrons in a bilayer system pierced by A = 1 (a) and A = 3 
(b) magnetic flux lines, leading to v = 2/1 and u — 2/3, respectively. 


tableau [A, A] is d\ = y^(A + 1)(A + 2) 2 (A + 3). For example, for A = 1 (N = 2) we have 
[1] © [1] = [2] © [1,1] (that is 4 x 4 = 10 + 6) so that [1,1] corresponds to the fully anti¬ 
symmetric irrep with dimension d\ = 6. In Refs. [28, .29] we have provided a “composite 
particle” picture (a term imported from the quantum Hall effect jargon [32]) to physi¬ 
cally interpret the [A, A] configurations as all possible ways of distributing 2A flux quanta 
among two identical fermions with spin 1/2 occupying two layers a (top) and b (bottom); 
see Figure [l] for an illustration. Here we reproduce the argument, which is the following. 
We have two electrons attached to A flux quanta each. The first electron can occupy any 
of the four isospin (layer, spin) states: |6t)> 6 j,), |aj/ and (aj,), at one Landau site of the 
lowest Landau level. Therefore, there are ( 4+ ^ _1 ) ways of distributing A quanta among 
these four states. Due to the Pauli exclusion principle, there are only three states left for 
the second electron and ( 3+ ^ _1 ) ways of distributing A quanta among these three states. 
However, some of the previous configurations must be identified since both electrons are 
indistinguishable and A pairs of quanta adopt ( 2+ ^ 4 ) equivalent configurations. In total, 
there are 



— (A + 3)(A + 2) 2 (A + 1) — d x 


( 11 ) 


ways to distribute 2A flux quanta among two identical electrons in four states, which turns 
out to coincide with the dimension d\ of the tableau [A, A]. In the BLQH jargon, this 
case corresponds to filling factor v = 2/A. The case v = 1/N (one electron) is related to 
the fully symmetric representation [N] with dimension (N + 1 )(N + 2 )(N + 3)/6, which 
coincides with the number of ways ( N // 3 ) of distributing N flux quanta among four isospiu 
states. 

The corresponding carrier Hilbert space di\ of the d^-dimensional unirrep of t/(4) is 
spanned by the set of orthonormal basis vectors 2j + m < A, \q a ,b\ < j} fulfilling 
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the resolution of the identity [28] 


A (A—m)/2 j 

I V \ j ’ m )( j ' m I 1121 

/ > / > / ^ I\<7a,gb I ’ A ' 

m=0 j = o ; l q a ,qb=-j 


where JT =0 .i means sum on j = 0, 1, §,... (the “angular momentum”). These are 

the U( 4) isospin-A analogue of the U( 2) spin-s orthonormal basis vectors |fc) of the pre¬ 
vious section (we are omitting the labels s and A from the basis vectors \k) and \q™ qf ), 
respectively, for the sake of brevity). The sixteen U( 4) generators r ^ are represented in 
T-L\ by sixteen operators T^ u whose matrix elements have been calculated in Ref. [25]. 
It is convenient to introduce the U( 2) 2 C U( 4) (antisymmetric Lorentz-like) generators 
by Mik = —\eikiToi and M 0i = — \T 3i , i,k,l = 1,2,3 and the isospin ladder oper¬ 
ators by T±n = (Tip ^ i^2/i)/2. We shall use the Einstein summation convention with 
Minkowski metric r/ /ti , = diag(l, —1, —1, —1) (unless otherwise stated) and denote, for 
example, T ±/i = T± = r] ,JlJ T± u . In the BLQH literature (see e.g. [27]) it is customary 
to denote the total spin S) = T 0i /2 = ie ik iM k i and pspin Pj = T i0 / 2, together with the 
remaining 9 isospin R k i = T k /2 operators. Angular momentum operators Su of layers 
£ = a, b (resp. pspin t = =p) are then given by 

Su = ^(TT 0i - T 3i ) = ^(M 0i t ^e m M kl ) = ^(^ - R i3 ), (13) 

so that S = Sb — S a and M 2 = = —4 (S 2 + S 2 ) = —2(S 2 + Rf) with R t = 

(Ru, R 2i , R 3i ), S = (Si, S 2 , S 3 ) and P = (Pi, P 2 , P 3 ). The linear and quadratic 7/(4) 
Casimir operators are Ci = T 0 o and 

C 2 = tw^r^-ir 0 2 0 = ir| 0 -lM 2 + i(f_^ + f + ^) 

= S 2 + P 2 + R 2 , (14) 


respectively, with R 2 
The basis vectors 

A lr*) = (2j+m 


= R\ + R% + R\ and 5^ the Kronecker delta. 

\ J q ’™ qb ) are eigenstates of P 3 , Su and M 2 with eigenvalues 

- S a \£J = q ,\= -8 j(j + 1)|£,). 


(15) 


We also have 


Ci 


9a,9b 


) = 2\\ j ’ m 


QaiQb 


), 


C 2 


1 j,m \ 

I 9a,9ft/ 


A(A + 4)1 


j,™ \ 

9a, Qb' 


(16) 


The generalization of spin-.s CS \z) in (|TJ) to isospin-A CS | Z) can be formally accomplished 
by, roughly speaking, substituting some scalars by 2 x 2 matrices. In particular, the 
complex scalar z E § 2 = U(2)/U(l) 2 is replaced by a 2 x 2 complex matrix Z E G 2 = 
U(A)/U(2) 2 , that is, Z = z^a^i with four complex coordinates z 1 * = tr(ZcqJ/2 e C. 
Sometimes we will also adopt the following decomposition for a matrix Z E G 2 


Z 



i- 






tan —e 1/3± , Vi 


f cos ^ — sin 

ysin ^e -1 ^ cos ^ J 


a,b, (17) 
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in terms of the eight angles 9 a> b,i!)± G [0,7r) and 4> a ,b,/3± G [0, 27t), where V a ^ represent 
rotations in layers 1 = a, b (note their “conjugated” character). 

Isospin-A CS can be obtained as an exponential action of ladder operators T± = 
on (namely) the lowest-weight state |g~=o^~= 0 ) and can be expanded in terms of 
the orthonormal basis vectors \i ,m n ) as 

'Qa,qb' 


I z) = 


e2 


itr(ZtT + )|0,0 


lo,o 


V^A yEA-m)/2 y^ 
Z— /m=0 Z-/o— q. A Z-/ 


j= o; 


J Qa,Qb 


. c/P’ m 
= —? ' Qa,qb 


(Z) \ j ’ m 


<la,qb' 


det(cr 0 + ZZt)V2 


det(cr 0 + ZZt)V2 


(18) 


with 




( y\ _ 


^ 2 J + 1 / A + 1 
A + 1 \2 j + m + 1 


A + 1 


m 


det (Zy‘VlJZ), 


(19) 


a set of homogeneous polynomials of degree 2 j + 2m in four complex variables z C gener¬ 
alizing the U( 2) monomials ipk(z) in ([!]). They are written in terms of the usual Wigner 
D- matrix m 


K,J Z ) 


(j + gJKj -&)! 

(j + ®>)Ki - ®>)! 


min(j+g a ,j-|-gi,) 

E 

fc=max( 0 , 5 a+ 9 t) 


(j+Qb\ ( j—Qb \ 
\ k J \k-q a -qb) 


„k j+q a —k j+q b —k k—q a —q b 
Z U Z 12 Z 21 Z 22 , 


for a general 2x2 complex matrix Z with entries z uv and angular momentum j. Homo¬ 
geneous polynomials (19) verify the closure relation (the Bergmann kernel K\) 


A (A—m)/2 j 

E E E (Z) = det( ff0 + Z'ZV = A' a (Z', Zt) 


m =0 j= 0 ;| qa,qb=-j 


( 20 ) 


from which we see that isospin-A CS are normalized but not orthogonal, with CS overlap 


(Z'\Z) 


Kx(Z',&) 

K x/2 (Z f ,Z^)K x/2 {Z, Z t)‘ 


( 21 ) 


Using orthogonality properties of the homogeneous polynomials tp 3 ^ (Z), a resolution of 
unity for isospin-A CS has been proved in [28], namely 1 = f G \Z)(Z\dfj,(Z, Z^), with 
integration measure 


d/x(Z, Z f ) 


i2 d x nZo d2z>x 

7 r 4 det(ao + Z^Z) A 


3 d x 
2 9 7T 4 


(cos'd.). 


cos d-) 2 df2 + df]-dQ a d^lb, 


( 22 ) 


where dC2± = sin ’d±dd±d(3± and dfle = sin Ocddidcpe (£ = a,b) are solid angle elements. 
This is the G2 generalization of the § 2 integration measure given after eq. (J5]) . 












3.1 A matrix generalization of the binomial distribution 

The probability of finding in a CS | Z) follows the multivariate distribution function 




i^irji 2 - 

2.7 + 1 




Qa,qb 




A + 1 


det(cr 0 + ZZ t) A 

A + 1 \ /A + 1\ det(ZZt)" 

2j + m + l) \ m ) det(cr 0 + ZZ^) X ' 


KJ 2 )?- (23) 


This is the Husimi function of a basis state | q™ qb ) (see later on Sec. 3.3). Let us study 
several particular cases. For example, using the parametrization (17), the Husimi function 

b ) = |oq) is given by the simple expression 


of the lowest-weight state y 


Qa,qbl 


^0,0/ y yf\ f 2^ + 

Qo,o( Z ’ Z ) = ( cos ~Y cos Yj ’ 


whereas for |^ ) = | 0 ’™) we have 

/A+i wA+n 

Q O ,m / y yf\ _ \m+l/ \ m / 

0,0 ^ i 


2 d + 2 
cos — COS 


~2~ 


A —m 


2 • 2 
sin — sm 


~2~ 


(24) 


(25) 


which has a binomial-like structure [similar to (|3])] with parameter cos 2 ^ cos 2 The 
general case (23) has a cumbersome explicit expression in terms of angles (17), but closed 
expressions can still be written for particular matrices Z. In fact, for diagonal Z — S = 

Co" £ ) ’ £± = tan % el/3± ’ and usin § that V q a ,qbY) = S qa,q b ¥+ qa Y qa we arrive at 


.2 0- 




,9b 


A +1 
A + 1 \2j + m + 1 

X-j-m 


2 $+ 2 1?_ 
X I cos — cos — 
2 2 


A + 1 
m 


2 d + 2 ^- XJ+m 

sm — sm — 

2 2 


tan 2 ^ N 
tan 2 


Other simplifications are also possible for the cases: 


V] j(Z) = (.:° + c 3 ) 2 '. P ; ; xZ) = (. z° - z 3 ) 2j , 
Vi-^Z) = (z 1 - i z 2 ) 2j , V j _ jtj {Z) = (z 1 + i z 2 ) 2j . 


(26) 


This matrix generalization of the binomial structure is also apparent from a Fock-space 
representation of isospin-A CS ( [18]) given in Ref. [28]. In this Fock-space picture, boson 
annihilation operators a and b in (|4j) are replaced by their matrix counterparts 


a = 





(27) 
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where 2 annihilate flux quanta attached to spin d and electrons 1 and 2 in layer t = a, 6 
in a BLQH context. If denote by a = /^"tr^a)^ and b = |?/ /it/ tr(cr At b)cr !y , then the CS 


| Z) in (18) can be written as a boson condensate of 2A flux quanta 


_i / d 2 E ±££ lV | o F> 

AI^XTI \ydet(ff 0 + ZF) / 


(28) 


and the Jordan-Schwinger boson realization of the sixteen 17(4) isospin operators r MI/ 
in terms of eight boson operators a and b is written as T IIV = tr (ZW^Z) with Z = 
(Z h Z 2 ) = (ab) a “two-fermion compound”. We have proved in [29] that wave func¬ 
tions are antisymmetric (fermionic character) under the interchange of the two electrons 
(Zi,Z 2 ) (Z 2 ,Zi) for A odd, and they are symmetric (bosonic character) for A even. 


3.2 Operator symbols: BLQH order parameters 

Now we comment on the analytic picture provided by the CS (Fock-Bargmann) represen¬ 
tation \P(Z) of any state |//i) € Hx, defined as 

*(Z) = K x/2 (Z,Z'){Z\il>). (29) 


For example, the basis states | ijj) = \ 3 j™ qh ) are represented by the homogeneous polynomi¬ 
als tpij™ (Z) of degree 2 j + 2 m in four complex variables z* 1 , whereas general states \ij)) 
are represented by polynomials 4 '/(Z) of degree < 2A. This makes the space ?ix(& 2 ) 
of polynomials \H(Z) a reproducing kernel Hilbert space with measure dfi\(Z, Z^) = 
dfjb(Z, Z^)/K\(Z, Z^). Inside this CS picture, isospin 1/(4) operators T /IU are represented 


by differential operators T, w . Writing Z = 


z 2 


= ZuZ 11 , 


= 


d/dz A 


and dfj, = 


d/dZ/j, = d M , these differential operators have the following expression [we use the nota¬ 
tion commented before (13) 


Ad 


T 30 = 2(^ M a M - A), — dp, T +Il = z 2 T-f, - Zn%oi 


(30) 


so that T^ v d>{Z) = Kx{Z,Z^){Z\T llv \ , ij}). As already commented in section [2J this ana¬ 
lytic representation reduces the computation of (technically elaborate in principle) matrix 
elements like (Z\Tfj U \Z') to simple derivatives 


(Z|T„„|Z'} 


7;,„AA(Z,Z't) 

/u(Z,Zt)/u(Z',Z't)’ 


(31) 


with K\(Z,Zd) — (1 + zJ J, 8 llu z lu + z tJ "7'i llu z' u ) x . Here we want to analyze the operator 
symbols (T^) = (Z\T^ U \Z) = K^(Z, Z Jf )T flu K\(Z, Z^), which turn out to be 


/rp \ _ r\ \ —l+det(X^^) 

( 730 / — ZA det ^o +z tz) 1 

Irr\ \ _ r) \ Z/z+det (Z^)Z/j, 

W-M/ — ZA det(o-°+ZtZ) ’ 


(M \ - 2A 

\ 1V1 ^1 ZA det(a°+ZtZ)’ 

I rp \ _ r) \ z /1 +det(Z)z fl 

\ J +/7 — ZA det (<T°+ZtZ) ■ 


(32) 
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We can check that (T +fl ) = det(Z)(T_ /i ) - ^(T 30 ) and = z^T-J) - z u {T_ IJ ). In 

the same way, we can also compute symbols of quadratic (and more general) operator 
combinations (Tj^T^v) = K\ 1 (Z, Z^)T i _ lv (T^ v 'K\(Z 1 Z^)). For example, U (2) 2 -invariant 
(i.e., commuting with M /w ) quadratic operators: T| 0 , M 2 and 

l( T+/tT M + T _ /jT m) = 4 (P 2 - P 2 2 - P? + Rl) 

i(T + ^-T_^) = 4(P,P 2 - R { • /f 2 ) (33) 

T_ m T^ = P 2 - P 3 2 + R 2 - R 2 3 - 4P 3 , 


have the following symbols 


(7f 0 ) 


A—1 /rp \ 2 | a \ 1+det Z) irp r T 1 f Jj \ _ 4A(A+1) det(Z^) /rp r T 1 L l \ _ 4A(A+1) det(^) 

~A~ \ 30/ ~r ^ A det(o-°+ZtZ)’ ~ det(a°+ZiZ) ’ \ i +^ 1 +/ — det(<j°+ZtZ) ’ 

(pr?) = 2 (a + \ h±M^021 _ ^l(r 3 „) 2 ) , (34) 

= ^(Tj 2 ,,) + yy(T 30 ) 2 - 2A(A + 2) 

= |<2?o> + 4(23o) + 2(pTf) - 2A(A + 4), 


where the last equality is a consequence of the quadratic Casimir (|14|) and the fact that 

[T +Ii ,f- V ] = 
(J34 ) only de- 

2 |M 2 - 


(T +/J T^) = (TL^T^) + 4(T 3 o), which comes from the commutation relation 
^T 30 +2M^ and the antisymmetry of M^ v . Note that all quadratic symbols 
pend on U( 2 ) 2 invariant quantities like tr(ZZ^) = |£ + | 2 + |£_| 2 and det (ZZ^) 
In [22] we have studied the interlayer imbalance (P 3 ) = (T 30 )/2 and its 
(P|) — (P 3 ) 2 in a CS. Here we want to analyze general pspin (P) and spin 
Before, we should stress that non-zero fluctuations a\ = {A 2 ) — {A) 2 of spin 
pspin (H = P) exist in a CS, so that (A 2 ) ^ [A) 2 in general. For example, it 


fluctuations 
(S) symbols. 
(H = S ) and 
is interesting 


to compare the Casimir (14) symbol 


(S 2 ) + (i 32 ) + (R 2 ) = A(A + 4), 


(35) 


with the quantity 

(S ) 2 + (P ) 2 + (R ) 2 = A 2 , 
so that a non-zero “variance sum”, written as: 


+ < 7)3 + — A(A + 4) — A 2 — 4A, 


(36) 

(37) 


exists and is proportional to A. Here we shall restrict ourselves to the study of ( S ) 2 and 
(P) 2 , leaving the analysis of fluctuations aside. 

The ground state of a BLQH system at one Landau site and at filling factor v = 2 
is a coherent state | Z). Depending on the strength of interlayer tunneling, bias voltage, 
Zeeman, Coulomb, etc, interactions, the BLQH system can appear in three phases: Spin 
or ferromagnetic (pspin singlet) phase, pspin (spin singlet) phase and canted or antifcr- 
romagnetic phase The three phases are characterized by the order parameters (S ') 2 
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Phase: 

Spin 

Pspin 

Canted 

Order 

parameter: 

(£) 2 = A 2 
(P) 2 = 0 

(S') 2 = 0 
(P) 2 = A 2 

(S ') 2 7 ^ 0 

<P) 2 ^ 0 


Table 1: Spin and pspin coherent mean values in the three BLQH phases. 


and (P) 2 with values in Table [ij In the literature (see e.g. [27]), only the case A = 1 is 
considered for this phase diagram study. Here we shall determine the regions of Cp (i.e., 
the values of Z) associated to each of the three phases: 


1. Spin phase: ( S ) 2 = A 2 . From eq. (36), we see that this condition is equivalent to 
(P) 2 = 0 = (R) 2 . Remember that (Si) = i end(Mki). Using the mean values in 


eq. (32) and the parametrization (17), we can write (S') 2 = A 2 /(d + , 9_)g(9 a , 9 b , A</>), 
with 


/($+,$_) = -(costP — cos$ + ) 2 , 
g(9 a ,9 b , A0) = -(1 — cos 9 a cos 9 b — cos A0 sin# a sin 9b), 


(38) 


and A <j) = <j) a — (j) b - The range of / and g is the interval [0,1]. The condition 
(S) 2 = A 2 then implies /($+,$_) = 1 = g(9 a ,9 b , A0). The condition /(i? + ,i?_) = 1 
is fulhllcd for ($_,i9 + ) = ( 7 r, 0) or (0, 7 r). The condition g(9 a , 9 b , A0) = 1 is fulhllcd 
for example for 9 a + 9 b = tt and A 0 = 7r, although some other possibilities exist. 

2. Pspin phase: (P ) 2 = A 2 . Remember that (P 3 ) = and 

= (T+q) + (T-q) A a[tr(Z)(l + det(Zt)] 


det(cro + ZZ^) 


= (7+q) - (T-q) A 3[tr(Z)(l - det(Zt)] 


2i 


det (cr 0 + ZZ^) 


(39) 


From eq. (36), we see that this condition is equivalent to (S') 2 = 0 = (R) 2 . In 
particular, from (38), we see that ( S ) 2 = 0 if / = 0 or g = 0. We have that 


/ = 0 if d_|_ = IP = d and that g — 0 when, for example, 9 a = 9 b and Acj) = 0 
(other solutions are possible). In these particular cases we have that (P ) 2 = A 2 (l — 
sin 2 [(/P - p + )/2] sin 2 7 ?), which fulfills (P ) 2 = A 2 when = 0 (that is, Z = 0 ) or 
when /3_ = /3 + (i.e., Z = tan(i?/2)cr 0 ). 

3. Canted phase: (S') 2 7 ^ 0 7 ^ (P) 2 - This condition is accomplished almost everywhere 
in G 2 , except for the zero-measure sets defining spin and pspin phases. 

A deeper study of BLQH Hamiltonians and their phase diagrams will be done else¬ 
where. 
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3.3 Husimi function and Wehrl’s entropy: a conjecture 

The Husimi function of a given state \ifj) G Hx is 


Qi,(z,z f ) = \(z\i >)\ 2 = 




det(a 0 + ZZt) A ' 


( 40 ) 


For unit norm states (^|^) = 1 , Q </, is normalized according to f G> Qp(Z, Z^)d/i{Z, Z^) = 
1 , with measure (22). At the end of Section [2] we have commented on Wehrl’s and Lieb’s 
conjectures on the entropy of ordinary and spin-s CS, respectively. In Ref. [31], Wehrl 
conjectured that any Glauber (harmonic oscillator) coherent state has a minimum Wehrl 
entropy of W min = 1. This conjecture was proved by Lieb [3B], who also conjectured that 
the extension of Wehrl’s definition of entropy to spin-s CS | z) will yield a minimum entropy 
of Wmin = 1 — l/(2s + 1) (note that W m \ n —> 1 for high spin values). This conjecture has 
been recently proved by Lieb and Solovej in [35], and an extension to SU(N) coherent 
states seems to be also plausible, at least for the fully symmetric representation (complex- 
projective C p N ~ l CSs). 

Here we extend the Wehrl’s definition of entropy of a state ^ to Grasmannian isospin-A 
CSs | Z) as 


W^ = 


' G 2 


Qip(Z, &)lnQ^Z,^)dp(Z, Zt). 


(41) 


This quantity is a measure of the area occupied by ip in phase space. We conjecture that 
this area is minimum when ^ is a CS (i.e., |?/>) = | Z') for any Z' G G 2 ) and the minimum 
value is given by 


Wrr, in = 


4A(2 + ^ + A)(2-^f + A) 
4(1 + A)(2 + A)(3 + A) 


= 4 


1 + A 2 + A 3 + A 


(42) 


It is direct to prove that this minimum value is attained in particular by the lowest-weight 
state \^>) = (go) (he., Z' = 0), whose Husimi function is given by the simple expression 

0,A 


(24), and by the highest-weight state | ijj) = | 0 ’ 0 ) (he., Z' = cx)), whose Husimi function i 


is 


given by (25) for m — A. For other particular CSs | Z'), we have numerically checked that 


they also attain the minimum Wehrl’s entropy (42). A formal proof of this conjecture will 
be investigated elsewhere. 

Remember that, for f/(2)/C(l) 2 -CSs we had RAm —> 1 for high spin s values. Now, 
for [/(4)/C(2) 2 -CSs we have H Anin —> 4 for high isospin A values. This limiting entropic 
value coincides with the complex dimension rfc(P) of the corresponding phase space P, 
namely: dc(§ 2 ) = 1 for the two-sphere and dc(& 2 ) = 4 for the Grassmannian. 

Another measure of the area occupied by i/j in phase space is the Husimi second 
moment: 

M^= [ Q 2 (Z,Zt)dp(Z,Zt), (43) 

J g 2 

which has to do with the so called “inverse participation ratio” (IPR). Broadly speaking, 
the IPR measures the spread of a state |-0) over a basis {|i)}f =1 . More precisely, if pi 
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is the probability of finding the (normalized) state \ip) in |i), then the IPR is defined as 
M,/, = If If/) only “participates” of a single state |£ 0 ), then p io = 1 and = 1 

(large IPR), whereas if \ifj) equally participates on all of them (equally distributed), p t = 
1/d, Vi, then M./ = 1/d (small IPR). Therefore, the IPR is a measure of the localization 
of | ij>) in the corresponding basis. For our case, the Husimi second moment (43) measures 
how close is "0 to a coherent state | Z'). We also conjecture that the maximum Husimi 
second moment 


= 


(2 + A) 2 (3 + A) 


4(1 + A)(l + 2A)(3 + 2A) 


1 

16 


1/2 45/32 3/32 


1 + A 1 + 2A 3 + 2A 


(44) 


is attained when if} is a CS. Remember that, for £7(2)/£7(l) 2 -CSs we had M max —> 1/2 
for high spin s values. Now, for £7(4)/£7(2) 2 -CSs we have M max —* 1/16 for high isospin 
A values. This limiting IPR value coincides with 2~ D , with D = d<c(P) the complex 
dimension of the corresponding phase space P. 


4 Conclusions and outlook 


We have analyzed the structure and properties of a previously introduced family of CS 
labeled by points Z ( 2x2 complex matrices) in the Grassmannian G 2 = £7(4)/£7(2) 2 . 
Applications to bilayer quantum Hall systems at fractions of filling factor v = 2 have 
been commented all along the paper in order to exemplify some abstract mathematical 
constructions. It is interesting to compare the Grassmannian CS in eqs. (18) and (28) 
with the Bloch CS in eqs. (JTJ) and Q, together with the CS overlaps (21) and 


respectively. We perceive a similar structure between them, although the Grassmannian 
case is more involved and constitutes a kind of “matrix Z generalization of the scalar z”. 
The matrix noncompact £7(2,2) version of £7(1,1) CS has also been explored by us in 
(351 05], and applications to massive conformal particles and pairing systems have been 
commented. 

We have also extended Wehrl’s definition of entropy to Grassmannian CS and con¬ 
jectured a lower bound for it (resp. an upper bound for the Husimi second moment). 
The study of localization properties of Hamiltonian eigenfunctions is also an important 
subject and we believe that our result will be relevant here. Actually, CS in general are 
essential to study the classical, thermodynamic of mean-field limit of algebraic quantum 
models undergoing a quantum phase transition, providing order parameters and the cor¬ 
responding phase diagram. As already stated at the introduction, the Husimi function 
and its entropic measures have proved to be useful in the characterization of quantum 
phase transitions of several molecular and optical models [3HH23> 2P,I2D 3D] and topolog¬ 
ical insulators m The use of Grassmannian CS as variational states to approximate the 
ground state energy of bilayer quantum Hall systems is in progress. We hope that they 
will provide new insight into the subject. 
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